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Abstract 

Double differential production cross sections of 7r~ and ir + mesons 
and the number of participating protons have been measured in central 
Au+Au collisions at 1 A-GeV. At low pion energies the tt~ yield is strongly 
enhanced over the n + yield. The energy dependence of the ir~ /tt + ratio 
is assigned to the Coulomb interaction of the charged pions with the pro- 
tons in the reaction zone. The deduced Coulomb potential increases with 
increasing pion cm. energy. This behavior indicates different freeze-out 
radii for different pion energies in the cm. frame. 

Relativistic heavy ion collisions are a unique tool to study nuclear matter at 
high densities and temperatures [Q. At bombarding energies around 1 A-GeV, 
baryonic densities of 2-3 times normal nuclear matter density (po = 0.17 fm -3 ) 
are expected to be reached during a time span of 10-15 fm/c (see e.g.: @). 
Thereafter, the high pressure built up in the reaction zone drives the system 
apart and particles freeze out according to their mean free path. A major 
experimental challenge is to determine the baryonic density in the different 
stages of the reaction. 
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This determination requires to measure both the number of baryons inside 
the reaction zone and its time-dependent volume. The number of participat- 
ing nucleons is derived from the projectile spectator charges ||, whereas the 
measurement of their spatial extent is more complex. Two-particle correlations 
(tttt,KK,pp) can be used to extract information on the size and lifetime of the 
source at freeze-out [Q, || . 

As an alternative approach, presented in this Letter, the Coulomb potential 
of the nuclear fireball is extracted by studying its interaction with oppositely 
charged pions. The strength of the Coulomb field depends on the number of 
participating protons and their geometrical configuration (resp. the freeze-out 
radius) at the time of pion emission. At a beam energy of 1 A-GeV pion produc- 
tion mainly proceeds via the excitation and decay of the A33(1232)-resonance. 
In this case the overall 7r~/7r + -ratio is determined by the iV/Z-ratio of the col- 
liding nuclei via isospin conservation. The variation of the 7r~/7r + -ratio as a 
function of the pion energy is attributed to the Coulomb field of the charge 
in the reaction volume which de- or accelerates the pions according to their 
charge. Previously, the influence of Coulomb effects on the emitted pions has 
been studied for different reaction systems and energies without deducing 
the fireball's Coulomb field. For example, in Ne+NaF collisions at 0.4 A-GeV 
a strong peak in the 7r~/7r + -ratio has been found at laboratory angles close to 
zero degrees which was explained by Coulomb interaction of the charged pions 
with the projectile spectator 

In this Letter we report on the measurement of double differential produc- 
tion cross sections of positively and negatively charged pions emitted in central 
197 Au+ 197 Au collisions at 1 A-GeV incident kinetic energy impinging on a target 
of 1.93 g/cm~ 2 . The experiment was performed with the Kaon spectrometer || 
at the heavy-ion synchrotron SIS at GSI. The spectrometer covers a momentum- 
dependent solid angle of AQ = (15-35) msr with a momentum resolution of 
dp/p ~ 0.01 over the full momentum range. The momentum resolution is domi- 
nated by multiple scattering. The momentum acceptance is Pmax/Pmin ~ 2 for 
a given magnetic field setting. The measured laboratory momenta vary between 
0.156 GeV/c and 1.5 GeV/c. The spectra were obtained at laboratory angles of 
(44 ± 4) degrees. This angular range corresponds to a range of normalized ra- 
pidity of 0.43 < y/ybeam < 0.74 with transverse momenta between 0.12 GeV/c 
and 1.2 GeV/c. The particle identification is performed by the reconstruction 
of the trajectory and by the determination of the particle velocity with two 
time-of-flight arrays. The collision centrality is determined by means of the hit 
multiplicity of charged particles in a 84-fold segmented scintillation detector 
close to the target covering polar angles between 12 and 48 degrees. The cor- 
responding number of participating protons is derived from the summed charge 
of projectile spectator fragments measured at forward angles by a 380-fold seg- 
mented hodoscope. The method is discussed in ref. || . In order to minimize the 
influence of spectator matter (nucleons not contributing to the reaction zone) 
we have selected collisions with the highest charged particle multiplicity and 
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pion emission around mid-rapidity. These central collisions which were selected 
represent (14 ± 4) % of the total reaction cross section. The corresponding 
number of participating protons is determined to be < Z par t >= 110 ± 8. 

Figure [l] shows the n~ and n + cross section d 2 a / '(dE*™ d^l c .m.) as a function 
of the kinetic energy E™£ in the center-of-momentum frame for central Au+Au 
reactions. The ir~ yield clearly exceeds the ir + yield at low kinetic energies 
whereas at higher energies these two yields are close to each other. Similar 
results are obtained by the FOPI collaboration, yet over a smaller momentum 
byte @. 

In order to determine the energy-integrated tt /7r + -ratio R eX p = 
(d<j(ir~)/dn)/(da(ir + )/d£l) we parameterized the pion spectra by the sum of 
two Maxwell-Boltzmann distributions and extrapolated the fit into the non- 
measured energy range (for details see ref . Jll[ ) . 

The overall 7r _ /7r + -ratio of R eX p — 1-94 ± 0.1 agrees well with the ratios 
derived from the isospin decomposition of individual nucleon-nucleon collisions 
at the given beam energy: Using the parameterization of ref. [ fl2| one obtains 
a value of Ri SO = 1.90. The corresponding value from the isobar model which 
assumes a pion production exclusively via a isospin 3/2 resonance is Ra = 1.95. 
The agreement between experimental and theoretical values indicates that the 
7r~/7r + -ratio reflects the N/Z asymmetry of the colliding system. 

Figure || shows the 7r~/7r + -ratio for Au+Au as a function of the pion kinetic 
energy in the cm. system. The 7r~/7r + -ratio decreases from 2.8 at low pion 
energies to a nearly constant value of 1.1 for pion energies above 0.4 GeV. 

For the following analysis we assume that the 7r _ /7r + -ratio is independent 
of the pion energy if Coulomb effects are disregarded. This assumption is sup- 
ported by transport model calculations (IQMD j fL3| , BUU iQ) which find - 
within the error bars - a constant 7r~/7r + -ratio of ~ 1.9 for all pion energies 
if the Coulomb interaction is switched off. The BUU calculation also includes 
baryonic resonances heavier than the A33(1232) resonance. 

In order to extract the strength of the Coulomb force we use a static ap- 
proximation for the Coulomb field. This assumption is justified as we consider 
only pion kinetic energies above 0.14 GeV. Their velocities (/3 > 0.87) are sig- 
nificantly larger than the expansion velocity of the nuclear fireball ((3 ks 0.4). 
(Recently, the influence of an expanding source has been studied for charged 
kaons showing a negligible effect at higher kinetic energies Jlq| .) According to 
ref. JTq| , the Coulomb force distorts the pion spectra by modifying both the 
kinetic energies of the particles and the available phase space. One derives the 
following relations for the measured cross sections at a given momentum p 



°(P) = = °o(po(p)) 



d 3 Po 



d A P 



(1) 



with the undisturbed momentum fio and the undisturbed differential cross sec- 
tion <7o(po(p)) taken at the momentum po(p). 
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The Jacobian in cq. (|l|) (Coulomb phase-space factor) is evaluated assum- 
ing spherical symmetry with the relativistic energy-momentum relation and by 
using the identity pdp — EdE: 



d 3 Po 



d A p 



p E dE 



pEdE 



As the pion spectra are not described by a single Maxwell-Boltzmann distribu- 
tion, the undisturbed differential cross section is parameterized by d 3 ao/dpQ oc 
exp(-0(E o ) ■ E Q ) with p- 1 = 32 MeV + 7.4 • 1CT 2 • E Q . This function has 
been obtained as the average of the functions (3(E ) obtained for the measured 
7r~ and 7r + spectra. For inclusive reactions the function f3 determined in this 
way describes well the curvature of neutral pion spectra measured in the same 
reaction Jl7| . 
With (c=l) 



E± = E ± Vcoui; P± = yjE% - m 2 (3) 
one writes for the spectra of charged pions 

<r±(p±) = f± ■ cq(Po(p±)) • P ° ° ■ (4) 

P±E± 

The plus (minus) sign refers to ir + (ir~), and f± reflects the weights due to the 
isospin coefficients. With = /_//+ the tt~ /tt + ratio can be written as 

_ a -{pi) _ tot a (p (pl)) p + E + 

' <y+(pl) tso ' Mpo(pI)) ' p-E- [) 

It is worth noticing that for p — > oo the value of R does not necessarily approach 
unity. 

In a first attempt we determine the Coulomb energy Vcoui using eq. ([s]). 
However, with the assumption of a constant Coulomb potential Vcoui one cannot 
describe the ir~ /n + ratio in the entire energy range. According to eq. (^|), a 
7r~/7r + ratio of around 1.1 as measured for pion kinetic energies around 0.6 GeV 
requires a Coulomb potential of 22 MeV. This value, however, is at variance 
with the 7r~ /tt + ratio at lower pion energies (see the solid line in fig. ^). Here, 
the observed ratio clearly exceeds the solid line suggesting a smaller value for 
the Coulomb potential at lower pion energies. Other attempts to describe the 
measured n~ / n + ratio with a constant Vcoui failed as well |l8|, . 

One can esti- 

mate the freeze-out radius of high-energy pions (E* ™ = 0.6 GeV) from the 
determined Coulomb potential of 22 MeV. We assume that the pions are emit- 
ted from the surface of a charged sphere, containing all participating protons. 
The assumption of surface emission is justified by the large pion re-absorption 
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probability inside the reaction volume. The Coulomb potential is then given by 
Vcoui = Zpart-e 2 - < r- 1 >. With < Z part >= 110±8 and with V Co ul = 22 MeV 
one obtains an effective radius of 

r eff =<r^ >~ 1 = (7.2 ±1.1) fm. (6) 

According to p e ff = iA par t/{A'Kr^^) with the number of participating nucleons 

A part this value corresponds to an effective density of p e // = (1-1 -0,3) ' Po- 

In the following step we determine the Coulomb potential for lower pion 
energies. The disagreement between the solid line (calculated with a constant 
Coulomb potential ) and the data in fig. ^ is a signature for a different Coulomb 
field acting on low-energy pions. In order to describe the measured ratio R 
over the entire energy range we vary in eq. (^) the Coulomb potential Vcoui 
as a function of the pion kinetic energy. The phase space factor is modified 
accordingly. The resulting dependence of Vcoui on is shown in fig. ||: Vcoui 
changes from about 22 MeV for high-energy pions to less than 10 MeV for pion 
energies below 0.2 GeV. At even lower pion kinetic energies pions have de Broglie 
wavelengths comparable to the size of the emitting source; this problem requires 
treatment in terms of quantum dynamics which is beyond the scope of this letter. 

A value of 10 MeV for pion energies below 0.2 GeV corresponds to an effective 
freeze-out radius of r e // = (16 ± 2) fm and an effective density of p e // = 
(0.1 ± 0.03) • po- The reduction of Vcoui therefore indicates a more dilute charge 
distribution at freeze-out for low-energy pions and consequently larger effective 
freeze-out radii. 

A dilute nucleon distribution is expected for the late stage of the reaction 
whereas higher densities are reached in an earlier stage. This interpretation 
is in agreement with microscopic transport models which predict that high- 
energy pions (transverse momenta pr > 0.5 GeV/c) are emitted during an early 
stage of the collision (t < 15 fm/c) where the central baryon density is above 
2 • po while low-energy pions are emitted later during expansion |f20| , |2l|] . This 
scenario is also supported by the experimental observation, that the multiplicity 
of high-energy pions > 0.45 GeV) increases more than linearly with the 

number of participating nucleons |2^| . This behavior is a signature of pion 
production via multiple baryon-baryon collisions which occur more frequently 
at higher densities. For the system Ar+KCl at an incident energy of 1.5 A-GeV 
a variation of the source size with the pion energy has been found previously 
from studying two-pion correlations [p3[ . 

Recently 7r _ /7r + -ratios have also been measured at much higher bombarding 
energies. In Au+Au collisions at 10.7 A-GeV the 7r~/7r + -ratio is found to be 
1.5 for low transverse momenta pr and it approaches unity at higher px [ p4[ . In 
Pb+Pb collisions at 158 A-GeV the tt~ /7r + -ratio exceeds one only for very small 
values oipT- This weak energy-dependence is explained by Coulomb interaction 
with a rather small co-moving charge of Z e j f = 40 [^5|, ^6) . Two features of the 
measurements at the higher incident energies differ from our observations: at 
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high bombarding energies the number of pions per baryon is significantly higher 
and due to charge conservation the overall 7r~/7r + -ratio tends to be close to 
unity [^7| . The small energy dependence reveals a rather weak Coulomb effect 
(low number of charges Z e ff) likely due to a fast disintegration of the source. 
In contrast, the large variation of ir~ /tt + from 2.8 to about 1.1 at high pion 
energies as found in the Au+Au system at 1 A-GeV (see fig. ||) indicates a 
strong Coulomb effect due to a slowly expanding source (compared to the pion 
velocity) with a large effective charge Z e f f . 

In summary, we have measured 7r~ and n + spectra in central Au+Au colli- 
sions at 1 A-GeV together with the number of participating protons. For this 
heavy system we found a strong enhancement of the ir~ yield over the ir + yield 
at low pion kinetic energies. Based on the assumption that the n~ /n + ratio 
is given by Coulomb interaction and isospin conservation, the Coulomb poten- 
tial of the reaction volume is determined. With the additional experimental 
information on the number of participating protons, effective source radii are 
derived. They are found to be smaller for high-energy pions than for low energy 
pions. This result can be interpreted as a consequence of a pion emitting source 
which expands: high-energy pions are emitted more likely in the dense and hot 
phase of the reaction while low-energy pions leave the system after expansion 
at a more dilute phase. 
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Figure 1: Production cross section of negatively and positively charged pions 
from central collisions of the reaction system 197 Au+ 197 Au at an incident beam 
energy of 1 A-GeV and at an emission angle of 6lab — (44 ± 4) degrees. 
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Figure 2: n~ /ir + ratio for the reaction system 197 Au+ 197 Au as a function of 
the pion kinetic energy for central collisions. The solid (dotted) curve shows 
the calculated ratio for a fixed Coulomb energy of 22 MeV (10 MeV). The error 
bars reflect statistical errors only; a systematic error of ±5% has to be added. 
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Figure 3: Variation of |Vcw| with the pion kinetic energy as deduced from the 
measured tt~ /n + ratio for central 197 Au+ 197 Au collisions. 
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